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Abstract 


The  time  and  spatial  evolution  of  the  cathode  region  of  a  planar 
discharge  in  helium  has  been  modeled  using  numerical  solutions  of  the  electron 
and  positive  ion  continuity  equations  and  Poisson's  equations.  Single  step 
ionization,  electron  and  ion  drift  and  diffusion,  and  electron-positive  ion 
recombination  are  included  using  rate  and  transport  coefficients  based  on 
experiment.  The  model  is  applied  to  high  pressure,  high  current  discharges 
such  as  found  in  electron-beam  or  photon  initiated  lasers.  The  cathode  fall 
is  found  to  approach  steady-state  voltage  drop  and  sheath  thickness  are 
compared  with  experiment  and  previous  theory  using  appropriate  scaling  laws. 
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I.  INTRODUCTION 


The  time  and  spatial  evolution  of  ionization  growth  in  planar  discharge  in 
helium  has  been  studied  using  electron  and  positive  ion  rate  equations, 
Poisson's  equation,  and  experimental  transport  and  rate  coefficients.  A 
numerical  method  is  described  for  tracing  theoretically  the  transient  behavior 
of  a  rapidly  developing  gaseous  discharge  where  the  space-charge-distorted 
field  is  very  important.  We  undertook  this  study  of  the  cathode  fall  region 
as  an  initial  step  toward  the  solution  of  the  growth  of  discharge 
instabilities  in  a  three-dimensional  geometry.  In  addition,  there  appear  to 

1  O 

be  only  very  limited  calculations1-  the  time-dependent  growth  of  this 
boundary  layer.  This  paper  presents  numerical  calculations  of  the  quantities 
characteristic  of  cathode  fall  phenomena  and  a  comparison  of  the  numerical 
results  with  experiments  and  other  theories.  Section  II  describes  the 
fundamental  theory  and  numerical  scheme  used  in  our  calculations.  The  results 
of  the  calculations  are  shown  in  Section  III  and  are  discussed  in  Section  IV. 
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II.  THEORY 


A.  Basic  Equations 

The  planar  discharges  considered  here  are  the  discharges  developed  on 
an  area  with  a  diameter  large  compared  with  the  gap  length.  They  can  be 
treated  theoretically  by  a  one-dimensional  model.  The  electron  and  positive 
ion  rate  equations,  Poisson's  equation,  and  experimental  transport  and  rate 
coefficients  are  used  for  tracing  the  growth  in  time  and  space  of  a  one¬ 
dimensional  helium  discharge  where  the  electric  field  is  distorted  by  the 
space  charge  of  the  electrons  and  ions.  The  continuity  equation  for  electrons 
is 


9ne 

■W  *  1  •  (Vl>  =  vne  -  °W  0) 

where  ng  and  n+  are  the  electron  and  positive  ion  densities,  t  is  the  time,  v 
is  the  ionization  coefficient,  and  a  is  the  recombination  coefficient.  The 
electron  velocity  u  can  be  expressed  in  terms  of  drift  and  diffusion  velocity 
as 


V(n  D) 

u  =  -pE - - -  (la) 

e 

where  E  is  the  electric  field,  D  and  y  are  the  diffusion  and  mobility 
coefficients  of  the  electrons,  respectively . 

The  continuity  equation  for  positive  ions  is 

3n+ 

lt“  +  7  *  (n+Ej  =  vne  '  anen+  *  W 
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The  ion  velocity  u,  can  be  expressed  as 


y+ = 


y+£. 


(2a) 


where  y+  is  the  mobility  of  positive  ions.  The  diffusion  of  the  ions  has  been 
neglected  because  of  their  heavy  mass.  In  these  calculations,  the  ion  drift 
flux  is  always  greater  than  the  ion  diffusion  flux  by  at  least  a  factor  of 
100. 

Poisson's  equation  describing  the  space-charge-distorted  fields  is 


3E 

32 


(3) 


where  e  is  the  electronic  charge,  is  the  permittivity  of  free  space  and  z 
is  the  position  in  the  axial  direction. 

The  coefficients  D,  v,  a,  u+  and  y  must  be  taken  from  the  experimental 
measurement  or  electron  transport  coefficient  calculations.  For  atomic  helium 
they  are  expressed  in  terms  of  the  ratio  of  electric  field  to  gas  density  E/N 
in  the  unit  of  V  -  crrr  and  shown  in  Eq.  (8)  of  Reference  4.  Some  of  them  have 
been  modified  because  of  the  larger  range  of  E/N  throughout  the  discharge  gap 
(i.e.  10~18  V-cm2  <  E/N  <  5  x  10~18  V-cm2).  Those  used  in  the  numerical 
calculations  are  listed  for  convenience  as  follows: 


yN  =  2.3  x  1022  (cm-V-sec)"1 

y+N  =  2.9  x  1020  (300/Tn)0,26  (cm-V-sec)'1 

a  =  [1.13  x  10"8  +  4.22  x  10"29  N  T^]  300/T(°k)  cm3  sec'1  (4) 
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v/N  =  4.6  x  10 5 f E/N)  (1  +  5  x  1014  E/N)  exp(-7  x  10“16  N/E)  cm3  sec"1 


D  =  ykT/e 


where 


T  =  300  +  3  x  1020  (E/N)/  |l  +  [5  x  1010  (E/N)2/3 ]  }  '  . 

Here  D  is  the  longitudinal  diffusion  coefficient,  k  is  the  Boltzmann  constant 
and  T^j  is  the  gas  temperature  in  Kelvin. 

An  allowance  has  been  made  for  the  time  required  for  electrons  emitted 
from  the  cathode  to  attain  the  ionization  rate  appropriate  to  the  local  E/N 
value  by  requiring  the  electrons  to  pass  through  a  potential  of  40  V  before 
ionization  begins.  This  approximation  is  an  extension  of  the  observations  for 
a  uniform  electric  field5  to  the  present  spatially  varying  electric  field.  In 
practice,  the  following  smooth  function  has  been  used  for  the  ionization  rate 
coefficient. 


/IV10 

U  %)  for  0  <  V  <  40 

|  1-  G^)10  for  40  <  V  <  80  . 

B.  Boundary  Conditions 

At  the  cathode,  secondary  electrons  are  produced  by  the  bombardment  of 
ions  and  photons 

re(0)  =  y+r+(0)  +  Yprp,  (5) 
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where  l\  represents  the  flux  of  various  particles  i.  y+  and  Yp  are  the 
secondary  coefficients  representing  the  probability  of  the  production  of 
secondary  electrons  at  the  cathode  per  incident  positive  ion  or  photon.  The 
effect  of  photons  is  to  enhance  the  production  of  secondary  electrons  and 
speed  up  the  formation  of  the  cathode  fall  at  the  early  times  in  discharge 
which  will  be  demonstrated  in  Section  III.B. 

At  the  anode,  we  have  zero  density  for  both  electrons  and  ions 

ne(d)  =  n+(d)  =  0  *  (6) 

where  d  is  the  distance  between  two  electrodes.  The  voltage  across  the 
discharge  gap  is  V  so  that  the  boundary  condition  for  Eq.  (3)  is 

d 

f  Edz  =  V.  (7) 

0 

Equation  (1)  is  a  second  order  partial  differential  equation  with  two 
boundary  conditions.  Equations  (2)  and  (3)  are  first  order  partial 
differential  equations  with  one  boundary  condition.  Therefore,  we  have  a 
complete  set  of  three  partial  differential  equations  that  can  be  solved 
numerically  for  the  three  unknown  ne,  n+  and  E.  Note  that  the  boundary 
condition  in  Eq.  (5)  is  a  combination  of  Neumann  and  Dirichlet  boundary 
condition  which  has  the  form  An  +  B  (3n/3z)  =  C  where  A,  B  and  C  are 
constant.  An  imaginary  point  outside  the  boundary  has  to  be  introduced  in  the 
finite  difference  frame.  For  ion  equation,  we  can  use  one-sided  difference 
approximation  for  the  cathode  boundary. 
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C.  Current  Continuity  and  External  Circuit 

The  large  potential  drop  between  the  cathode  and  the  negative  glow  is 
normally  called  the  cathode  fall.  It  is  one  of  the  most  important  quantities 
describing  the  glow  discharge.  In  order  to  obtain  relations  between  the 
cathode  fall  Vc,  current  density  J  and  thickness  of  the  cathode  fall  dc  and 
compare  the  calculated  quantities  to  the  experimental  measurements,  an 
external  circuit  with  a  constant  external  resistance  £2  and  a  constant  voltage 
source  Vg  have  been  added  in  the  discharge  problem.  The  effects  of  circuit 
inductance  and  capacitance  has  been  neglected  so  that  we  have  the  simple 
circuit  equation 


V0  =  V  +  JoA, 


(8) 


where  J  represents  the  total  ionization  current  density  which  can  be  - 
calculated  by  the  following  energy  balance  method. 

The  power  fed  by  the  voltage  supply  V  across  the  discharge  gap  is  eqaual 
to  the  sum  of  the  power  delivered  to  the  gas  and  the  rate  of  change  in 
electrostatic  energy,  i.e., 

d  _  d  enE2 

Jv  -  /  I  •  J  <iz  +  -ft  /  -f-  dz  •  (9) 

0  0  c 

Equation  (9)  can  be  further  simplified  for  constant  V  by  using  Green's  theorem 


d 

Jd  =  e  /  (ngu  +  n+u+)dz  .  (10) 
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It  will  be  noticed  that  Eq.  (10)  relates  the  energy  input  to  the  gap  to 

the  integral  of  the  energy  dissipated  per  unit  volume  and  include  changes  in 

the  energy  stored  in  the  electric  field,  i.e.,  the  displacement  current.  In 

the  numerical  calculations,  for  each  time  step  we  use  Eq.(10)  to  calculate  J 

first  and  then  use  Eq.  (8)  to  adjust  V  to  insure  consistency  of  V  with  the 

2 

assumed  Vg  and  oA.  The  value  of  1  flm  for  oA  has  been  used  in  the 
calculations  throughout  the  paper  exclusively. 


D.  Numerical  Scheme 

A  Crank-Micholson  finite  difference  scheme'7  which  is  implicit  both  in  time 
and  space  has  been  used  and  a  predictor-corrector  step  has  been  added  to  take 
into  account  the  nonlinear  coefficients  in  the  equations.  Use  of  nonuniform 
spatial  step-sizes  has  been  found  very  effective  in  studying  the  time 
development  of  the  boundary  layer  in  the  cathode  fall  region.  Equations  (1) 
to  (3)  can  be  generalized  by  the  following  equation 


an 

at 


S  n  +  + 


€ 


(ID 


where  6,  p,  y  and  £  may  be  functions  of  z  and  n. 

Using  Taylor  series  expansion  or  Lagrange  interpolation  formula,  the 
first  and  second  spatial  derivatives  can  be  approximated  by  a  centered 
difference  in  space  as 


dn 

dz 


Vj-1 +  Vj  +  Cinj+1 
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and 


8 


d2„ 


=  A„n.  ,  +  B0n.  +  C„n. 


dz 


2  2  j-1  2  j  2  j+1 


The  Crank-Nichol son  scheme  yields  for  the  general  Eq.  (11) 


nt+l  t 
n .  -n . 


nt+ 1  t 
n .  +  n . 


t+1.  t 


t+1.  t 


j  J  _  „t  ,  "j  "j  ,  tf  n"ri+  n  ,  tP  n  +  n" 

~E - 5j  +  6j  - 2 -  +  ^j  5x  - 2 -  +  Yj  5xx - 2 -  * 


(12) 


where  we  use  the  operator  notation 


6  nt  =  A,  n1:  .  +  8^*  +  C,n* 
x  1  J-1  1  J  1  J+1 


5xxn  ~  A2nj-1  +  B2nj  +C2nj+1 


Equation  (12)  along  with  boundary  conditions  [Eqs.  (5)  to  (7)]  can  be  solved 
for  each  unknown  by  inverting  a  tridiagonal  matrix. 

In  order  to  have  the  scheme  completely  implicit  and  stable,  the 
coefficients  8,  n,  y,  5  have  to  be  centered  in  time.  We  then  use  a  predictor- 
corrector  iterative  scheme  with 


nj  "  nj  _  -v+1/2  pv+1/2  nj>  +  nj  -vfl/2  „  n V+1  +  n* 

E  C  +  8  - 2 -  +  71  5x - 2 - 


+  y^1/2S 

XX  2 


and 
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We  could  iterate  to  convergence  by  choosing  the  number  of  iterations  p  at  each 
time  step  so  that 

AP-1 

Max  n.  -  n.  <  s  . 

j  J  J 

t+1  *p 

For  suitable  e,  we  could  stop  the  iteration  and  set  n.  =  n..  Usually  a 

J  J 

stable  scheme  takes  no  more  than  two  iterations  for  each  time  step. 

As  discussed  in  Appendix  A  the  diffusion  term  in  Eq.  (1)  is  very  important 
in  the  calculation  in  conjunction  with  the  numerical  instabilities.  The 
expression  for  D  in  Eq.  (4)  has  been  chosen  carefully  for  the  best  fit  to  the 
experimental  observations.  Experience  suggests  that  in  order  to  solve  the  set 
of  Equations  (1)  to  (3)  without  numerical  instability  we  have  to  modify 
Eq.  (la)  as 

Vn 

U  =  -  yE  -  D  (13) 

ne 

where  the  diffusion  coefficient  D  has  been  outside  of  the  gradient.  Eq.  (13) 
has  been  used  throughout  the  paper  exclusively.  Furthermore,  in  the  time 
evolution  of  the  cathode  fall,  the  electron  equation  becomes  very  stiff  in  the 
early  time  around  0.07  ysec  for  the  case  run  without  introducing  the  external 
circuit.  Because  of  this,  the  set  of  partial  differential  equations  [Eqs.  (1) 
to  (3)]  has  been  solved  with  the  diffusion  term  [second  order  term  in  Eq.  (1)1 
increased  by  a  factor  of  2.  The  enhanced  diffusion  term  is  used  to  damp7  the 
high  frequency  terms  in  the  solution  without  affecting  the  desired  solution 
too  much  (see  Appendix  A). 
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III.  NUMERICAL  RESULTS 


A..  Initial  Conditions 

A  stable  difference  scheme  is  not  unduly  sensitive  to  small 
perturbations  so  that  an  initial  condition  can  be  chosen  arbitrarily.  As  in  a 
double  discharge  laser  the  discharges  considered  here  are  initiated  by  an 
electron  and  ion  pair  whose  density  varies  smoothly  along  z  and  has  a  peak 
value  of  3  x  10^  cm-3.  This  initial  ionization  with  Lorentz  profile  as  shown 
in  Figure  1  has  the  following  simple  relation  with  z 

n(z,0)  =  3  x  1010  [1  +  (z/d  -  0.5)2]"16  cm-3  .  (14) 

An  alternate  initial  condition  is 

n(z,0)  =  1.2  x  10^ (z/d)  (1  -  z/d)  cm-*3.  (15) 

Where  the  same  value  of  3  x  10^®  cm-3  was  chosen  at  the  center  of  the  gap  in 
Eqs .  (14)  and  (15).  Calculations  are  shown  for  a  gap  length  of  2  cm,  a  helium 
density  N  of  3  x  1019  cm-3,  temperature  TN  of  300°K,  and  initial  E/N  of  1.8  x 
10"  V-cm  .  A  y+  of  0.09  has  been  chosen  and  the  effect  of  photons  has  been 
neglected  in  the  calculations  shown  except  in  Figure  7.  It  is  known  that  the 
quantities  characteristic  of  cathode  fall  phenomena  do  not  depend  on  the 
length  of  the  discharge.  The  effect  of  gas  heating  is  neglected  at  this  stage 
so  that  Tjy  and  N  are  assumed  independent  of  time  and  space. 

Figure  2  shows  that  the  electric  field  is  reduced  by  charge  separation  in 
the  region  of  maximum  initial  ionization  at  the  very  earlay  time  (e.g., 

0.02  usee).  The  density  difference  between  electrons  and  ions  is  very  small 
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in  the  plasma  region.  The  field  near  cathode  increases  as  the  electrons  are 
swept  toward  the  anode.  The  significant  ionization  occurs  due  to  the  increase 
of  electric  field  in  the  cathode  sheath.  This  large  source  of  ionization 
causes  a  remarkable  increase  in  the  positive  ion  density  while  the  electrons 
are  moving  away  from  the  cathode  sheath  continuously.  The  field  increases 
near  the  cathode  until  it  becomes  high  enough  to  cause  an  ionization  wave  to 
move  toward  the  cathode.^  This  is  illustrated  in  Figure  3  for  two  different 
times.  Note  that  the  figures  show  only  the  small  region  within  0.03  cm  of  the 
cathode.  The  electric  field  near  cathode  becomes  high  and  forms  the  cathode 
fall.  At  t  =  0.12  usee,  as  shown  in  Figure  4,  the  current  density  is  3.4 

O  O 

amperes/cnr  and  the  thickness  of  cathode  fall  is  about  2.4  x  10  cm.  It  is 
noted  that  Figures  2  to  4  show  results  obtained  using  initial  profile  (Eq.  (14)) 
and  boundary  condition  (Eq.  (7))  with  constant  V.  As  mentioned  early  the  diffusion 
coefficient  D  has  been  enhanced  by  a  factor  of  2.  The  sequence  of  electric 
field  and  charge  distributions  in  the  cathode  fall  region  is  essentially 
independent  of  gap  length  with  the  time  scale  decreasing  from  ~200  ns  for  a  3 
mm  gap  and  100%  overvoltage  to  ~60  ns  for  a  7  cm  gap  and  50%  overvoltage. 

In  order  to  compare  our  results  with  other  theories,  the  circuit  equation 
(8)  has  been  added  in  the  calculations  from  this  point  on.  Figure  5  shows  the 
transient  behavior  of  the  discharges  with  nA  =  1  ohm  -  m2  for  initial  condition 
(Eq.  (15)).  The  similar  results  for  initial  condition  (Eq.  (14))  except  in  the 
very  early  time  suggests  that  we  can  choose  any  arbitrary  initial 
conditions.  The  distributions  of  electric  field,  electron  and  ion  densities 
change  drastically  prior  to  0.1  nsec  due  to  the  propagation  of  the  ionization 
wave.  It  then  takes  a  rather  longer  time  to  reach  the  steady-state  solution 
as  shown  at  t  =  0.4  usee.  We  see  little  difference  between  the  curves  at 
0.3  ysec  and  at  t  =  0.4  ysec.  Following  convention,  the  cathode  fall  distance 
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and  cathode  fall  voltage  are  determined  by  extending  the  straight  line  portion 
of  the  electric  field  to  the  point  where  the  electric  field  is  zero.  This  is 
illustrated  in  Figure  6  which  shows  the  values  of  Je,  J.,  E,  and  J  of  Figure  5 
at  t  =  0.4  ysec  as  a  function  of  the  distance  from  the  cathode.  Note  that 
J  =  J0  +  is  the  total  current  density  and  is  independent  of  the  distance 
from  the  cathode. 

At  the  cathode  there  are  few  electrons  and  the  current  is  carried  mainly 
by  ions;  in  the  negative  glow  region  where  the  field  is  very  small,  the  ion 
and  electron  numbers  are  more  equal  but  because  of  their  greater  mobility  the 
electrons  carry  the  current.  Experimentally  field  measurements  are  used  to 
define  dc  as  the  point  where  the  extrapolated  field  becomes  zero9  and  then 
determine  Vc.  While  none  of  the  theories  has  yet  been  able  to  account  for  the 
experimentally  observed  linear  change  of  field10  in  the  cathode  fall  region, 
the  calculations  presented  here  have  shown  the  linear  decrease  in  field  from 
cathode  to  the  edge  of  the  glow. 

B.  Effect  of  Electron  Reflection  of  Cathode 

In  this  subsection  we  examine  the  sensitivity  of  the  solution  to  the 
boundary  condition  of  the  cathode.  The  Milne  boundary  condition11  with  the 
inclusion  of  cathode  injection  and  electric  field  can  be  expressed  as 

ve  %  dnP 

ne  T  {1  "  r)  "  Yn+V+E  =  D  dz“  ‘  VE  »  (16) 

where  the  electron  mean  velocity  v0  =  (2  kT/m)1//2,  r  is  the  electron 
reflection  coefficient  and  y  is  the  secondary  coefficient  from  cathode.  The 
left-hand  side  of  Eq.  (16)  is  the  net  flux  of  electrons  arriving  at  the 
cathode  from  an  imaginary  boundary  one  mean-free-path  into  the  gas,  while  the 
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right-hand  side  is  the  same  flux  as  given  by  the  usual  transport  coefficient 
formul ation . 

We  choose  r  =  0.2  and  y  =  0.3  and  do  the  similar  calculations  as  in  A 
with  new  boundary  conditions.  The  results  shown  in  Figure  7  are  very  similar 
to  Figure  5.  Two  different  boundary  conditions  do  not  give  any  significant 
difference  in  the  time  evolution  of  the  cathode  fall  region.  The  small 
difference  for  the  electric  field  in  the  cathode  sheath  is  due  to  the 
difference  in  the  secondary  coefficient  y.  Note  that  y  is  defined  differently 
in  both  cases.  Fig.  8  shows  the  growth  of  the  electron  and  ion  current 
density  as  well  as  the  voltage  drop  through  the  gap.  We  can  see  the  current 
has  been  saturated  and  reached  the  steady  state  for  the  case  with  external 
circuit.  The  difference  of  current  density  and  voltage  drop  in  the  early  time 
is  due  to  a  different  choice  of  initial  condition.  A  140  Volt  cathode  fall 
was  found  at  t  =  0.4  ysec. 


C.  Role  of  Metastables  and  Resonance  Radiation 

The  rate  equations  for  helium  triplet  (23S)  metastables  T  and  singlet 
(2*S)  metastables  S  are-*-2 


3T 

3t 


=  k  Tn  N 
eT  e 


°T  32T 


3Z 


«T3n2t 


(17) 


and 


3S 

3t 


kesne  N  +  N  2  “  kS2NS  "  kS3^ 

9Z 


(18) 


where  Ker  is  the  rate  coefficient  for  excitation  of  species  y  at  unit  gas 
density,  ky3  is  three-body  collision  frequency  for  conversion  species  y  into 
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metastable  molecules,  is  two-body  colllslonal  frequency  for  species  y,  Dy 
is  the  diffusion  coefficient  at  unit  gas  density  for  species  y.  The  current 
density  of  metastable  atoms  reaching  the  cathode  by  diffusion  can  be  neglected 
on  our  short  time  scale.  Also,  since  the  metastable  molecules  do  not  radiate 
they  do  not  produce  a  photocurrent . 

Now  the  K$2  term  in  Eq.  (18)  is  collision  induced  radiation  so  that  the 
current  of  photons  arriving  at  the  cathode  is 


rs  = 


kS2 

2 


S  dz 


(19) 


Note  that  rg  includes  the  fraction  (1/2)  of  nonresonance  photons  produced  in 
the  gap  which  are  directed  toward  the  cathode. 

The  differential  equation  governing  the  density  of  resonance  radiation 
can  be  written  as  follows1^ 


3t  ~  keRneN  "  kRlN  R  "  kR2N  R  ~  kR0R 
d 

-  AR  +  A  f  R(z' ,t)K( |z-z' | )dz '  ,  (20) 

0 

where  A  is  radiative  transition  probability  for  the  resonance  transition,  keR 
and  Krq  are  frequencies  of  excitation  at  unit  gas  density,  KR1  and  KR2  are 
two-body  and  three-body  collision  frequency  respectively ,  K(jz-z'|)  is 
transmission  function  for  resonance  radiation.  The  integral  of  K(|z-z'|) 
gives  the  probability  that  the  resonance  photons  emitted  at  z  will  be  absorbed 
at  a  point  to  the  left  of  z  =  0  in  the  absence  of  the  cathode  or  the 
probability  the  resonance  photons  will  travel  as  far  as  the  cathode  before 
absorption. 
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If  we  use  the  Biberman  approximation. 


12 


then 


d 

R  -  /  R(z‘ ,t)K( J z-2 ' | )dz'  =  R[G(z)  +  G(d-z)]. 
0 


?  1/4 

We  use  the  approximation  G(z)  =  l/[  16  +  31  ( -nk „z)  j  ,  where  k  =  3tt/ X 
and  Xg  is  the  wavelength  of  the  resonance  line.  The  current  of  photons 
rR  becomes 


d  0 

rR  =  A  /  R(z,t)dz  f  K( | z-z 1 | )dz‘ 
0  -<» 

=  A  /  R(z,t)dz  f  K(y)dy 
0  z 

d 

=  A  f  R(z,t)G(z)dz  . 

0 


8oth  r  and  rn  contribute  to  r  in  Eq.  (5).  The  coefficients  used 
s  R  p  ^ 

from  Table  II  and  Fig.  3  of  Reference  12  are  listed  below. 


k  =  2.8  x  105  E  (1 - E/N  exp(-3.7  x  1016  N/E)  sec-1 

es  9  x  10'iO 

k  „  =  6  x  10"15  sec"1 

S2 

-32  -1 

k  .  <  10  sec  1 
s3 

KeR  =  1.38  x  105  E  exp  (-3.8  x  10"16  N/E)  sec"1 

-12  3  -1 

kR^  =  2  x  10  crri  sec 


(21) 

come 
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kR2  =  9  x  10"32  cm6  sec"1 
kRg  =  1.9  x  10^  sec"1 
Xq  =  5.84  x  10"6  cm 

A  =  1.8  x  109  sec  .  (22) 

We  next  perform  calculations  similar  to  those  in  section  III. A,  i.e., 
with  initial  condition  (Eq.  (14))  and  boundary  condition  (Eq.  (5)).  A  yp  of 
0.05  has  been  chosen.  The  results  are  shown  in  Figure  9  and  should  be  compared 
to  Figure  5  in  order  to  see  the  effect  of  photons  on  the  evolution  for  the  cathode 
fall.  As  pointed  out  previously,  the  time  scale  for  the  formation  of  the  cathode 
fall  has  been  shortened  because  of  increased  production  of  secondary  electrons  at 
the  cathode  at  early  times.  Other  than  that,  the  effect  of  photons  can  hardly 
be  seen  in  the  later  stage  of  the  discharge  because  the  ion  term  becomes 
dominant.  A  cathode  fall  of  183  Volt  was  found  at  0.4  visec  for  the  conditions 
of  both  Figure  5  and  Figure  9. 
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IV.  DISCUSSION 


Although  the  mean  free  path  for  electrons  with  an  average  energy  of 
~5  eV  is  small  compared  to  the  cathode  sheath,  a  number  of  authors  have 
questioned  the  applicability  of  the  continuum  model  implicit  in  Eqs.  (1)  and 
(2).  As  a  partial  answer  to  this  question  we  calculate  the  energy  relaxation 
distance^  for  conditions  corresponding  to  those  of  Figure  4.  The  energy 
relaxation  distance  for  the  plasma  region  is  about  4  x  10~®  cm  which  is  very 
small  compared  to  the  dimension  of  the  plasma  region.  A  relaxation  distance 
of  3  x  10-5  cm  is  found  at  the  midpoint  of  the  cathode  fall  region.  Since  the 
electric  field  varies  by  only  4%  over  this  distance,  we  believe  that  the 
continuity  equations  are  reasonably  accurate. 

Theoretically  and  experimental ly  numerous  studies  of  the  cathode  fall 
region  have  been  made  to  explain  the  mechanism  of  the  cathode  phenomena  of  the 
glow  discharge.  Various  theories  have  been  derived  from  a  selection  of  the 
aforementioned  fundamental  equations,  i.e.,  equations  fl)  to  (3);  most  people 
use  Townsend's  a  for  the  rate  of  ionization  and  make  different  assumption 
about  u+  and  a  as  a  function  of  local  electric  field.  For  example.  Von  Engel 
and  Steenbeck10  omitted  continuity  equation  and  used  instead  the 
experimental ly  determined  distribution  of  the  field  in  the  cathode  fall 
region.  Similar  calculations  have  been  made  by  Ward-*-  except  that  he  had  to 
use  different  form  for  a  and  u+  in  order  to  get  good  agreement  with 
experimental  data.  Experimentally ,  VQ  has  been  defined  and  measured  in  a 
number  of  ways.  The  definition  of  Vc  depends  strictly  upon  the  method  used  to 
measure  it.  Each  method  has  certain  inherent  errors.  Many  experimenters  have 
obtained  results  not  agreeing  with  each  other,  and  often  not  consistent  and 
reproducible  in  their  own  apparatus.  However,  extensive  cathode  fall 
measurements  have  been  made  by  Guntherschul ze  and  co-workers.^  In  order  to 


18 


eliminate  the  heating  effect,  they  used  in  his  measurements  a  very  large 
massive  cathode  and  low  gas  pressures.^  Warren18  has  developed  a  refined 
electron  beam  technique  to  measure  the  field  in  the  cathode-fall  region.  They 
all  show  an  almost  linear  decrease  of  field  to  zero  with  distance  from  the 
cathode.  No  negative  fields  were  found  in  the  negative  glow.  He  further 
pointed  out  that  there  are  two  theories  of  the  cathode  region  which  give 
results  not  agreeing  well  with  experiment.  One  which  assumes  that  ions  moves 
across  the  cathode  fall  region  freely  (fee  fall  theory);  the  other  which 
assumes  that  ion  has  a  drift  velocity  u+  proportional  to  the  square  root  of 
electric  field.  All  of  their  results  for  the  cathode  fall  versus  current 
density  have  been  shown  in  Figure  10  where  we  have  compared  our  results  for 
the  steady  state  for  cathode  fall.17  We  see  that  the  calculated  cathode  fall 
voltage  drop  is  close  to  steady-state  values  measured  at  much  lower  gas 
densities.  Note  that  the  cathode  fall  is  a  strong  function  of  secondary 
coefficient  y. 

Essentially  the  factor  which  determines  Vc  is  y  which  depends  on  the  size 
and  energy  of  the  ions  and  photons  hitting  the  cathode,  on  the  cathode 
material  and  the  state  of  its  surface,  and  on  the  nature  and  pressure  of  the 
gas.  The  calculations  have  been  made  with  an  external  resistance  and  fixed 
voltage  supply  so  as  to  simulate  experiments  more  accurately.  Other  processes 
which  should  be  included  in  the  model  include  the  production,  destruction  and 
ionization  of  metastable  atoms  and  molecules  and  the  heating  and  motion  of  the 
neutral  gas.18  Obviously,  further  work  is  necessary  to  clarify  the 
relationships  among  the  various  theoretical  approaches  to  the  cathode  fall. 

At  present,  both  the  continuum  model  and  the  nearly  free  fall  approaches  seem 
to  be  able  to  use  their  methods  to  show  their  validity. 
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ELECTRICAL  FIELD,  volt/cm 


DISTANCE  FROM  CATHODE,  cm 


FIGURE  1.  INITIAL  DENSITIES  OF  ELECTRON  AND  ION  AND  ELECTRIC  FIELD  IN 
PARALLEL  PLANE  GAP 
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ELECTRON  AND  ION  DENSITIES,  cm 


DISTANCE  FROM  CATHODE,  cm 


FIGURE  2.  CALCULATED  ELECTRON  AND  ION  DENSITIES  AND  ELECTRIC  FIELD  AT 
t  =  .02  usee 
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ELECTRON 


ELECTRICAL  FIELD,  volt/cm 


0  0.01  0.02  0.03 

DISTANCE  FROM  CATHODE,  cm 

FIGURE  3.  PROPAGATION  OF  THE  IONIZATION  WAVE:  CALCULATED  CHARGE  DENSITIES 
AND  ELECTRIC  FIELD  AT  t  =  .07  ysec  AND  t  =  .08  ysec 


ELECTRON  AND  ION  DENSITIES,  cm 


ELECTRICAL  FIELD,  volt/cm 


FIGURE  4. 


CALCULATED  ELECTRON  AND  ION  DENSITIES  AND  ELECTRIC  FIELD 
AT  t  =  .13  ysec.  THE  CURRENT  DENSITY  IS  10  AMPERE/CM2 


NEAR  CATHODE 
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FIGURE  5.  THE  EVOLUTION  OF  THE  CATHODE  FALL  REGION  WITH  INITIAL  CONDITION 
(EQ.  (15))  AND  fiA  =  1  OHM-M2 


ELECTRON  (- - -)  AND  ION  ( - )  DENSITIES,  cm 


E(V/cm) 


FIGURE  6.  THE  VALUES  OF  Je,  J-j,  E,  and  J  OF  FIGURE  5  AT  t  =  0.4  nsec 
VERSUS  DISTANCE  FROM  THE  CATHODE 
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(zwo/v)f 


ELECTRICAL  FIELD,  ( - ) volt/cm 


8  x  104 


DISTANCE  FROM  CATHODE,  cm 


FIGURE  7.  CALCULATIONS  SIMILAR  TO  FIGURE  5  EXCEPT  USING  THE  NEW  BOUNDARY 
CONDITION  (EQ.  (16))  WITH  r  =  0.2  and  Y  =  0.3 
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ELECTRON  (- — )  AND  ION  ( - )  DENSITIES,  cm 


ELECTRICAL  FIELD,  (- - )  volt/cm 


DISTANCE  FROM  CATHODE,  cm 


FIGURE  9.  CALCULATIONS  SIMILAR  TO  FIGURE  5  EXCEPT  THAT  THE  EFFECT  OF  PHOTONS 
HAS  BEEN  CONSIDERED 
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FIGURE  10.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT  FOR  THE  STEADY  STATE  OF 
THE  CATHODE  FALL  VOLTAGE  VERSUS  CURRENT  DENSITY 
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18.  The  calculated  rate  of  heat  input  near  the  center  of  the  cathode  fall 
using  the  results  shown  in  Fig.  4  is  1  mW/cm  or  2.3  x  10  eV/ atom-sec. 
This  corresponds  to  a  temperature  rise  of  about  250  K  in  100  nsec.  A 
sound  wave  would  travel  through  the  cathode  fall  region  in  about  10  nsec. 
The  expanding  low  density  region  produced  by  heating  in  the  cathode  fall 
region  has  been  observed  experimentally .  For  example,  Crawford  and 
Phelps  (unpublished)  have  made  interferometric  observations  of  this 
region  and  Rogoff  has  calculated  some  of  its  properties. 
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APPENDIX  A 


STABILIZATION  OF  FINITE  DIFFERENCE  SCHEME 
BY  MEANS  OF  A  DIFFUSION  TERM 


First  consider  the  unstable  difference  scheme  for  the  hyperbolic  equation 


Ut  +  cUx  =  0 


Un+1  =  un  _  2  (un  -  un  ) 

j  j  2  [Uj+i  j-r 


f  Al) 


where  c  is  a  constant  and  X  =  cAt/Ax. 

The  Lax-Wendroff  scheme  for  the  same  equation  is 


(Jn 

J 


1  (Un 

2  (  j+1 


n 

j+1 


2Uj +  u3-i}  • 


(A2) 


It  is  stable  provided  Jxj  <1.  We  could  regard  this  scheme  as  an  obvious 

2  2 

difference  approximation  to  the  equation  IL  +  clJY  =  dJ  where  e  =  x  Ax~/2At. 
In  the  limit  as  the  mesh  spacing  approaches  zero,  our  differential  equation 
becomes  Uj.  +  cUx  =  0.  Therefore  we  have  added  a  diffusion  type  term  in  (A2) 
to  an  unstable  scheme  (Al)  and  thereby  stabilized  the  scheme.  Also,  we  might 
expect  the  solution  of  the  difference  scheme  to  converge  to  the  solution  of 
the  hyperbolic  equation. 

Second,  let  us  consider  the  Taylor  series 
Uj"+1  =  U(xrtn  ♦  at) 


At 


■  Uf*rtn>  +  ““t'VV  +i% 


utt<xj'tn>  + 


out3) 


(«) 


O 

For  the  hyperbolic  equation  =  -cUx  we  have  =  c  Uxx,  thus  (A3)  becomes 


Ujn+1  =  U;-CAtUx+^Uxx  +  0(At3)  . 


(A4) 


A-2 


2 

It  is  clear  that  we  even  improved  the  accuracy  from  0( At)  to  0 ( At  )  which 
gives  from  (Al)  to  (A2).  This  improved  accuracy  comes  from  the  Taylor  series 
expansion  which  produced  the  stable  difference  scheme. 

The  diffusion  term  regarded  as  a  filter,  the  removal  of  energy  in  the 
high  frequency  modes  can  be  done  to  pass  low  frequencies  and  block  high 
frequencies.  If  we  consider  the  equation 


* 


U.  +  cl)  =  aU 
t  x  xx 


Using  a  Fourier  series  the  solution  is 

U(x.t)  =  y  a.  e-a’2k2t  ei'k(x-ct)  . 
k=-»  K 

2  2 

If  the  term  air  k  t  is  quite  small,  then  the  solution  approximates  that  of  the 
hyperbolic  equation  U.  +  cU  =  0,  namely 

T#  X 

U(x,t)  =  y  a,  elirkfx*ct)  . 
k=-«  K 
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